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Abstract 

We analyze the potential of the CERN Large Hadron Collider (LHC) to study the structure 
of quartic vector-boson interactions through the pair production of electroweak gauge bosons via 
weak boson fusion qq — ► qqWW. In order to study these couplings we have performed a partonic 
level calculation of all processes pp — ► jje^n^vv and pp — ► jje^fi^uu at the LHC using the exact 
matrix elements at 0(a® m ) and 0(a^ m a 2 ) as well as a full simulation of the ti plus to 2 jets 
backgrounds. A complete calculation of the scattering amplitudes is necessary not only for a correct 
description of the process but also to preserve all correlations between the final state particles which 
can be used to enhance the signal. Our analyses indicate that the LHC can improve by more than 
one order of magnitude the bounds arising at present from indirect measurements. 
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I. INTRODUCTION 



Within the framework of the Standard Model (SM), the structure of the trilinear and 
quartic vector boson couplings is completely determined by the SU (2) L x U{1)y gauge sym- 
metry. The study of these interactions can either lead to an additional confirmation of the 
model or give some hint on the existence of new phenomena at a higher scale The triple 
gauge-boson couplings were probed at the LEP 0, 0] and are still under scrutiny at the 
Tevatron jj| through the production of vector boson pairs. However, we have only started 
to study directly the quartic gauge-boson couplings |3| • If any deviation from the SM pre- 
dictions is observed, independent tests of the triple and quartic gauge-boson couplings can 
give important information on the type of new physics responsible for the departures from 
the SM. For example, the exchange of heavy bosons can generate a tree level contribution 
to four gauge-boson couplings while its effect in the triple-gauge vertex would only appear 
at one-loop level, and consequently be suppressed with respect to the quartic one 0. 

At present the scarce experimental information on quartic anomalous couplings arises 
from the processes e + e~ — > W + W~ , y : Z77, ZZ^y, and W77 at LEP 3]. Due to phase space 
limitations, the best sensitivity is attainable for couplings involving photons which should 
appear in the final state. Photonic quartic anomalous couplings can also affect 77Z and 
77W productions at Tevatron 0, Q] and they will be further tested at LHC |s] and in the 
long term at the next generation e + e~ collider 0, 10 , [ill 3 Q • 

Purely electroweak quartic couplings W + W~W + W~ and W + W~ ZZ have not been di- 
rectly tested so far but will be within reach at LHC jl3 . HR It^ . In this work we 
study the potential of the LHC to probe them by performing a detailed analysis of the most 
sensitive channels that are the production via weak boson fusion (WBF) of W + W~ pairs 
accompanied by jets, i.e., 

p + p — > jjW + W~ — > jje^/j^uu , (1) 

and the WBF production of a pair of jets plus W ± W ± 

P + P^ jjW ± W ± jje ± fi ± uu . (2) 

We have only considered final state with different flavor leptons (e and /i) in order to 
avoid backgrounds coming from Z, 7 — > e + e~ or fi + [i~. The advantage of WBF, where the 
scattered final-state quarks receive significant transverse momentum and are observed in 
the detector as far- forward/backward jets, is the strong reduction of QCD backgrounds due 
to the kinematic configuration of the colored part of the event. 

There are previous studies of the quartic gauge boson couplings at the LHC. The earlier 
works relied upon the equivalence theorem [1^ or/and the effective H^-boson 

approximation . In Ref. 0| the full tree level calculation of the processes pp — > VV+ 
2 jets, with V = W ± ,Z° was presented. Here, we improve over these earlier works by 
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computing the full matrix element for all processes with the six fermion final states in ((TJ 
and (J2J) at 0(a® m ) and 0(a^ m a^). This includes the contribution from the resonant gauge 
boson pair production considered in Ref. jl?} as well as all the non-resonant contributions 
and their interference. We have also performed a full simulation of the ti background and 
evaluated the ti plus 1 and 2 jets backgrounds using the narrow width approximation for 
the top. 

The interactions responsible for the electroweak symmetry breaking play an important 
role in the gauge-boson scattering at high energies as they are an essential ingredient to avoid 
unitarity violation in the scattering amplitudes of massive vector bosons at the TeV scale 
[2fl| . There are two possible forms of electroweak symmetry breaking which lead to different 
solutions to the unitarity problem: (a) there is a particle lighter than 1 TeV, the Higgs 
boson in the standard model, or (6) such particle is absent and the longitudinal components 
of the W and Z bosons become strongly interacting at high energies. In the latter case, the 
symmetry breaking occurs due to the nonzero vacuum expectation value of some composite 
operators which are related with new underlying physics. 

We parameterize in a model independent form the possible deviations of the SM predic- 
tions for the W + W~W + W~ and W + W~ ZZ quartic gauge couplings in these two differ- 
ent scenarios as described in Sec. |H] In the first case we assume the existence of a light 
Higgs boson and consequently we are lead to dimension eight effective operators where the 
SU(2) L x U(l)y gauge invariance is realized linearly. We also contemplate the scenario 
where no new heavy resonance has been observed that leads to the gauge symmetry being 
realized nonlinearly by using the chiral Lagrangian approach. 

Valuable information on the possibility of new physics effects can also be gathered from 
the low energy data and the results of the Z physics; see Ref. j^] for a recent review. 
In particular they can constrain the possible deviations of the quartic gauge boson self- 
interactions from the SM predictions through their contributions to the electroweak radiative 
corrections 22]. For completeness we present in Sec. IHII the updated bounds on these effects 
from the global electroweak fit. 

Sections|TV|and|V|contain the details of the strategies proposed to reduce the backgrounds 
to acceptable levels while keeping the signal from the quartic gauge vertex. We find that the 
complete calculation of the scattering amplitudes is necessary to preserve all correlations 
between the final state particles which can be used to enhance the signal. We also study the 
precision with which the background rate in the search region can be predicted which is the 
ultimately limiting factor. 

Our final quantitative results on the attainable sensitivity at LHC are presented in Sec. I VII 
We find that LHC can improve by more than one order of magnitude the bounds arising at 
present from indirect measurements and it is able to test deviations with the size expected 
in the scenario in which no light Higgs boson is found and the gauge symmetry is realized 
nonlinearly. 



3 



II. THEORETICAL FRAMEWORK 



In this work we focus on the study of the structure of the weak quartic couplings con- 
taining H^'s and/or Z's. For the sake of simplicity we will consider effective interactions 
that do not contain derivatives of the gauge fields. With this requirement there are only two 
possible Lorentz invariant structures contributing to each of the four gauge boson vertices 

0™ w = g ap g^ [W+WjW+Wf] , Of w = g^g* [W+W+WjWf] , 
a wz = ga p g7 s [w+ZpW-Zs] , Of z = g^g^ 5 [W+WpZ,Z s ] , (3) 

O zz = Of z = O zz = 4*9* [Z a Z p Z^Z 5 ] , 

and the Lagrangian for the four gauge boson vertex will be 

c vvvv s c vv> vv> + c vv Q vv . (4) 

In the SM, SU(2) l gauge invariance and renormalizability imply that 

n ww _ ww _ 2 jvz _ _ 2 r wz _ n 2 r zz _ n /c\ 

c 0,SM _ C 1,SM _ 2 C 0,SM — 2 C 1,SM — H C SM — u \° ) 

C W C W 

where cw is the cosine of the weak mixing angle and g is the SU(2)l coupling constant. 

Conversely, if the SM is thought of only as an effective low energy theory valid up to 
the scale A, one expects deviations from Eq. (JSJ) even if we still retain the gauge symmetry 
group, the fermionic spectrum, and the pattern of spontaneous symmetry breaking (EWSB) 
as valid ingredients to describe Nature at energies E < A. In this case one can still write 
the Lagrangian for the four gauge boson interactions as Eq. (@J) but now the coefficients, Cq 
and Ci will be in general independent, and we can write 

c r=cs+p 2 A C r. (6) 

In the language of effective Lagrangians the deviations Aq will be generated by higher 
dimension operators parameterizing the low energy effect of the new physics. The order on 
the expansion at which these deviations are expected to appear depends on whether the low 
energy spectrum still contains a light SM-like Higgs boson responsible of EWSB or, on the 
contrary, EWSB is due to a heavy (or not fundamental) Higgs boson. 



A. Effective Operators with Linear Realization of the SU(2)l x U(l) Gauge Sym- 
metry 

We first assume that the low energy spectrum contains a light Higgs boson. In this case 
we chose a linear realization of the symmetry breaking in the form of the conventional Higgs 
doublet field In the usual effective Lagrangian language, at low energy we describe the 
effects of the new physics - which will manifest itself directly only at scales above A - by 
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including higher-dimension operators in the Lagrangian. The basic blocks for constructing 
the operators which can modify the four gauge boson electroweak vertices are the Higgs 
field, its covariant derivative -D M $, the SU(2) L field strength W^ u , and U(l)y field strength 



B^ v . The lowest order operators which can be built are of dimension six [23(. However 
dimension six operators which modify the four gauge boson vertices, affect either the two or 
three gauge boson couplings as well. Consequently they are better searched for, and severely 
constrained at present, by looking into those effects. 

The lowest dimension operators that modify the quartic boson interactions but do not 
exhibit two or three weak gauge boson vertices are dimension 8. The counting is straight 
forward: one can get a weak boson field either from the covariant derivative of $ or from 
the field strength tensor. In either case the vector field is either accompanied by a vacuum 
expectation value (VEV) of the Higgs field (v ) or a derivative. Therefore genuine quartic 
vertices are of dimension 8 or higher. There are only two independent dimension 8 operators 
without derivatives of the gauge fields (for further details see appendix |XJ) 
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When the Higgs field $ is replaced by its VEV, (J7J) and (JHJ) generate four gauge boson 
interactions as Eqs. and (jBJ) with 
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B. Effective Operators with Non-Linear Realization of the SU(2)l x £7(1) Gauge 
Symmetry 

If the electroweak symmetry breaking is due to a heavy (strongly interacting) Higgs 
boson, which can be effectively removed from the physical low-energy spectrum, or to no 
fundamental Higgs scalar at all, one is led to consider the most general effective Lagrangian 
which employs a nonlinear representation of the spontaneously broken SU(2)l £g> U(l)y 
gauge symmetry j^]. The resulting chiral Lagrangian is a non-renormalizable non-linear a 
model coupled in a gauge-invariant way to the Yang-Mills theory. This model independent 
approach incorporates by construction the low-energy theorems |2j| , that predict the general 
behavior of Goldstone boson amplitudes irrespective of the details of the symmetry breaking 
mechanism. Notwithstanding, unitarity implies that this low-energy effective theory should 
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be valid up to some energy scale smaller than Airv ~ 3 TeV, where new physics would come 
into play. 

To specify the effective Lagrangian one must first fix the symmetry breaking pattern. We 
consider that the system presents a global SU(2)l <£> SU(2)r symmetry that is broken to 
SU(2)c- With this choice, the building block 1 of the chiral Lagrangian is the dimensionless 
unimodular matrix field E(x), which transforms under SU(2)i ® SU(2)r as (2,2): 

= exp i , (10) 



where the tp a fields are the would-be Goldstone fields and r a (a = 1, 2, 3) are the Pauli 
matrices. The SU{2)l ® U(l)y covariant derivative of S is defined as 

= ^E + ^E-^SyB,. (11) 

Quartic vector boson interactions are generated at second order (p A ) in the derivative 
expansion j^. For simplicity we will consider only interactions which respect the custodial 
SU (2) symmetry. At this order, there are only two such operators usually denoted as 

4 4) = a 4 [Tr(^K)] 2 , (12) 
4 4) = a 5 [Tr(V,Vn] 2 , (13) 

where we defined = (D^H) EL These effective operators generate four gauge boson 
interactions as Eqs. (jlj) and (0) with 

Ac^ w = g 2 a 4: ee Ac 0i „o-iin Acf w = g 2 a 5 = Ac 1>no _ lin 



n z 9 2 _ A c o,no-iin wz _ g 2 _ Aci 

,no— lin i -, a \ 

AC - - -p AC 1 - WjT a 5 - — Z-j (14) 



a zz 9 , , Ac ,no— lin 4" Aci no — Hn 

Ac = - 3 -(a 4 + « 5 ) = -j 

III. LOW ENERGY CONSTRAINTS 

Valuable information on the possibility of new physics effects can also be gathered from 
electroweak precision data, measured mainly at the Z-peak by LEP1 experiments, but also 
including the W and top masses and other measurements. These data can be used to 
constrain the possible deviations of the quartic gauge boson self-interactions from the SM 
predictions as they contribute to the gauge boson self-energies at the one-loop level j^. 

Standard Model electroweak radiative corrections as well as universal new physics ef- 
fects enter in the predictions of these electroweak precision observables in three different 



We follow the notation of Ref. 



24]. 
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combinations usually named 61,62,83 26] (or S, T, and U j^), so in general 

£ i = £i,SM + £i,new ■ (15) 

Technically the procedure to obtain the contribution from the operators (J7j), (jHJ), (fT2|) . 
and (fT3*j) to the e's is the following: first we evaluate their contribution to the self-energies 
using dimensional regularization. Then, we keep only the leading non-analytic contributions 
- that is, the terms proportional to log(/i 2 ) - dropping all others. These contributions are 
easily obtained by the substitution 



2 A 
log 



2 



4 - d ~° M| ' 

where A is the energy scale which characterizes the appearance of new physics. 

With this procedure we found in Ref. 22] that for the operators (J7j), (jSJ), (f*T2|) . and (fT3|) . 

^2,new = ^3,new = and that only £i in ew is non-vanishing: 

^new= -^^(l + ^flog^, (16) 

£i,Bew- -^2^(1 + Ch/ )^- log ^, (17) 

where Ac, for the case of linear [non-linear] realization of the gauge symmetry are defined 
in Eq. © [Eq. d] 

Recent global analysis of the low energy and LEP data j^] yields 

6 X = (5.0 ± 1.1) x 1(T 3 , 

while the SM prediction is a function of m t , rrih, u s and a em . We use m t = 178 GeV, 
a s {M z ) = 0.119 and a em {M z ) = 1/128.88. 

For the case with a light Higgs boson of rrih = 120 GeV and a new physics scale A = 2 
TeV we find that at 99% CL 

- 5.2 < f x 10~ 3 < 9.0 , 
-13 < fi x 10~ 3 < 22 . (18) 

In models without a light Higgs boson, the gauge-boson contribution to S\ is infinite as 
a consequence of the absence of the elementary Higgs. On the other hand, we must also 
include the tree level effect due to the 0(p 4 ) operator which violates custodial SU(2) and 
which absorbs this infinity through the renormalization of the corresponding coefficient. If 
the renormalization condition is imposed at a scale A, we are left with the contribution 
due to the running from the scale A to Mz- Therefore, the SM contribution without the 
Higgs boson will be the same as that of the SM with an elementary Higgs boson, with the 
substitution ln(M#) — > ln(A). For A = 2 TeV we get the following 99% CL bounds 

- 0.32 < ct 4 < 0.085 , 

-0.81 < a 5 < 0.21 . (19) 
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IV. CALCULATION TOOLS 



We concentrate on the study of the structure of quartic vector-boson interactions through 
the production of W ± W Zf and W ± W ± in WBF, with subsequent decay to e/x pairs and 
neutrinos. The signal is thus characterized by two quark jets, which typically enter in the 
forward and backward regions of the detector and are widely separated in pseudorapidity, 
by a significant transverse momentum imbalance, and by a pair e ± /i T or e ± /i ± . 

Significant irreducible backgrounds can arise from QCD and electroweak (EW) processes 
which lead to the same final state 

V + V -> jje ± /i T w , jje ± /j, ± vv , 

where the jets arise from a gluon or light quark production. They include "resonant" pro- 
cesses with the production and subsequent leptonic decay of W ± W :f or W ± W ± pairs (on- 
or off-shell) accompanied by jets, and "non-resonant" processes containing only one or no 
W's in the s-channel. Furthermore for different sign final leptons, a large QCD background 
is expected from the production and subsequent decay of top quark pairs together with 0-2 
jets. 

The six-particle amplitudes for the signal and irreducible backgrounds are simulated at 
the parton level with full tree level matrix elements. The SM amplitudes are generated using 
Madgraph in the framework of Helas routines. The anomalous contributions arising 
from the effective interactions © and (j!4|) are implemented as subroutines and included 
accordingly. We consistently took into account the effect of all interferences between the 
anomalous and the SM amplitudes, and did not use the narrow-width approximation for 
the vector boson propagators. For the treatment of the finite-width effects in massive 
vector-boson propagators we use a modified version of the complex mass scheme (3(J in 
which we globally replace vector-boson masses m v with m v — imyTy without changing the 
real value of sin 2 W (3]], l^l • This procedure respects electromagnetic gauge invariance. 
We have also performed a full simulation of the ti background and evaluated the ti plus 1 
and 2 jets backgrounds using the narrow width approximation for the top quark. We took 
the electroweak parameters a em = 1/128.93, mz = 91.189 GeV, mw = 80.419 GeV, and 
?7itop = 174.3 GeV. The weak mixing angle was obtained imposing the tree level relation 
cos^vv = m w/ m z, which leads to sin 2 9w — 0.222. In our calculations we used CTEQ5L 
parton distribution functions 33^ 

The general expression for the total cross sections for the processes considered can be 
written as 

<J = crbck + g 2 (&c )a + g 2 {A Cl )at + g 4 (Ac ) 2 a 00 + g\Ac x ) 2 a n + ^ 4 (Ac )(A Cl )aoi , (20) 

where Aq for the case of linear [non-linear] realization of the gauge symmetry are defined 
in Eq. Q [Eq. (JHJ)]. o"b c k contains the contributions from all the backgrounds described 
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above while o"o and o\ contain the interference between SM and anomalous amplitudes. 
For the case of a linear realization of the gauge symmetry they contain the contribution of 
the light Higgs boson exchange, which is absent in the non-linear case. In either scenario 
the anomalous contributions a , a\, <Joo, <?oi, and an, as well as the EW contribution to 
Obck in the absence of a light Higgs boson, do not respect the unitarity of the partial-wave 
amplitudes (af) at large subprocess center-of-mass energies M ww |34(. For higher WW 
invariant masses, rescattering effects are important to unitarize the amplitudes. Taking 
into account this fact, we conservatively we impose in these cases the cut Mww < 1-25 
TeV, which guarantees that the unitarity constraints are always satisfied. This requirement 
corresponds to a sharp-cutoff unitarization j^j]. 

An important feature of the WBF signal is the absence of color exchange between the 
final state quarks, which leads to a depletion of gluon emission in the region between the 
two tagging jets. Thus one can enhance the signal to background ratio by vetoing additional 
soft jet activity in the central region [3^. Certainly, a central jet veto is ineffective against 
the EW backgrounds which possess the same color structure as the signal. For the QCD 
backgrounds, however, there is color exchange in the t-channel and consequently a more 
abundant production of soft jets, with p T > 20 GeV, in the central region [37[. The proba- 
bility of an event to survive such a central jet veto has been analyzed for various processes 
in Ref. 0], from which we take the veto survival probabilities 0.8 (0.3) for electroweak 
(QCD) processes. Moreover, at the high-luminosity run of the LHC there will be more than 
one interaction per bunch crossing, consequently there is a probability of detecting an extra 
jet in the gap region due to pile-up. In Ref. 39] it was estimated that due to pile-up the 
jet-veto efficiency for a threshold cut of px = 20 GeV is 0.75. Taking into account these two 
effects we obtain that the veto survival probabilities are 

= q g x o 75 = Q 6 ; pQCD = Q 3 x 75 = 225 . (21) 

Constraining quartic gauge boson couplings in the WBF processes pp — > jjefiuu is es- 
sentially a counting experiment since there is no resonance in the WW invariant mass 
distribution. The sensitivity of the search is thus determined by the precision with which 
the background rate in the search region can be predicted. In order to access the size 
of these uncertainties we have employed four different choices of the renormalization and 
factorization scales which we denote by: 



CI ^ F = ^ R = Upl n+v l n )/2- 



C2 fi° R = y (Prji +PTj 2 )/2 an d f^f = where s is the squared parton center-of-mass 
energy; 

C3 fM R = [M° F = y/PTjl PTj2] 

C4 a 2 s {ix R ) = a s (p Tjl ) a s (p T j 2 ) and fjP F = min(p Tjl , p Tj2 ). 
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Finally, we simulate experimental resolutions by smearing the energies (but not direc- 
tions) of all final state partons with a Gaussian error given by AE/E = 0.5 /y/E © 0.02 if 
\r]j\ < 3 and AE/E = \j\[E © 0.07 if |^| > 3 (E in GeV), while for charged leptons we 
used a resolution AE/E = 0.1/y/E © 0.01. We considered the jet tagging efficiency to be 
0.75 x 0.75 = 0.56 while the lepton detection efficiency is taken to be 0.9 x 0.9 = 0.81. 



V. SIGNAL AND BACKGROUND PROPERTIES 



A. Basic cuts 



We initially impose the following jet acceptance cuts 

p{ > 20 GeV , \7]j\ < 4.9 , (22) 

in order to have a well defined tagging jets. We also demand lepton acceptance and isolation 
cuts 



M < 2.5 , r] J m - m <r] t < r/ max 
ARtj > 0.4 , AR U > 0.4 (23) 



mill 



Pt > Pt 

where 77 min ( max ) is the minimum (maximum) rapidity of the tagging jets and p™ m = 100 (30) 
GeV for opposite (equal) charge leptons. Since the signal events contains undetectable neu- 
trinos that carry some transverse energy from the event, we also require a missing transverse 
momentum 

Pmissing > 30 GeV . (24) 

The tagging jets are usually well separated in rapidity in the signal, therefore we demand 
the existence of a rapidity gap between them 

\Vji ~ vM > 3.8 , rjji ■ r)j 2 < . (25) 



B. Additional cuts for pp — > jje fi T vv 



The production of opposite sign leptons exhibits a very large background due to the 
production of tt pairs in association with 0, 1, 2 jets. In the it process the 6-quarks produced 
in the t decays are identified as the tagging jets. We denote by tij and tijj backgrounds 
those events where the additional jet(s) is (are) identified as the tagging jet(s) while one of 
the 6-jets in tij and both in tijj are soft and central, tij and tijj events where one (or two) 
of the additional jet are not identified as the tagging ones contribute to the QCD radiation 
of the corresponding it and tij background and their effect is included in the gap survival 
probabilities |4p| . 
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FIG. 1: Normalized distribution of the transverse momentum of the charge leptons for ti (red solid 
line), tij (blue dashed line), tijj (black dot-dashed line), SM irreducible production (magenta 
dotted line), and anomalous W + W~ contribution <too (blue solid line marked "ano"). We assumed 
m h = 120 GeV and applied cuts (J2U)-(J2SJ) but for with a relaxed cut -p l T > 20 GeV. 



The relevance of the tighter cut on the transverse lepton momentum to suppress the 
different backgrounds in pp — > jje^^vv is illustrated in Fig. ^ In order to further reduce 
these backgrounds we make use of the fact that QCD processes typically occur at smaller 
invariant masses of tagging jets compared to EW processes. This is illustrated in Fig. El 
where we show the normalized invariant mass distribution of the tagging jets for the different 
backgrounds and the anomalous contribution a 00 for — > jje^ pTvv. Consequently, in order 
to further suppress the backgrounds we also require a large invariant mass of the tagging 
jets 

Mjj > 1000 GeV , (26) 

which mainly reduces the ti events but still leaves a large background from tij and tijj 
production. These events can be very efficiently suppressed by vetoing additional soft jet 
activity in the central region. Consequently, we impose that the event does not contain 
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FIG. 2: Normalized jet-jet invariant mass distribution for tt (red solid line), ttj (blue dashed line), 
tijj (black dot-dashed line), SM irreducible production (magenta dotted line), and anomalous 
VF + VF~ contribution o"oo (blue solid line marked "ano"). We assumed rrih = 120 GeV and applied 
cuts dHjHESJ). 

additional jets with transverse momentum larger than 20 GeV in between the tagging ones, 

4 < 20 GeV if r/ min < Vj < V U . (27) 



Additionally we notice that the azimuthal angular distribution of the charged leptons 
relative to each other in the SM is different than in the anomalous contributions. The e ± /i =F 
pairs from the decay of the W pairs produced via the effective interactions and (fT3j) are 
preferentially emitted in opposite direction from each other. This is shown in Fig. El where 
we plot the normalized distribution of the azimuthal angle between the electron and the 
muon. Thus we impose also the additional cut 

<f efl > 2.25 rd. (28) 
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FIG. 3: Normalized distribution of the e/i azimuthal angle difference for tij (dashed line), irre- 
ducible background (dotted line), and anomalous W + W~ production (solid line). We assumed 
m h = 120 GeV and applied cuts ((22M2Z1)- 

Finally we make use of the fact that the anomalous contributions arising from the effective 
interactions (jHJ) and (|14jl lead to a growth of the cross section for large WW invariant masses; 
see Fig. Consequently we define the signal region 

M™ w > 800 GeV , (29) 

where we define the transverse invariant mass as 

M? W = (^(p T y + < + y/ft + <) 2 - (PT + M 2 (30) 

where p'x is the missing transverse momentum vector, p^f is the transverse momentum of 
the pair e-fi and m efJi is the e/i invariant mass. 

In Table HJ we illustrate the effect of the above cuts for pp —>■ jje^^vv. In the lines 
marked IREDH — we take into account the full scattering amplitude for the irreducible 
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FIG. 4: Normalized distribution of Mrp for ttj (dashed line), irreducible background (dotted 
line), and anomalous W + W~ production (solid line). We assumed nth = 120 GeV and applied 

cuts gZ HZfl . 

backgrounds. We separate the electroweak and QCD part of these backgrounds in order to 
show the effect of the veto survival probabilities. As illustration of the signal loss due to 
the imposed cuts, we also include the cross section for the anomalous term o"oo- From this 
table, we can see that the largest background is the tt + n jets production by three order 
of magnitude when we apply only the acceptance and tagging cuts. However, after cuts the 
dominant backgrounds are ttj and EW irreducible processes. 

C. Cuts for pp — > jje^/j. ,vv 

In the first column in Table HP we give the cross sections for pp —>■ jje ± /i ± z/^ after the 
basic cuts fl2"2"j) - (}2l)j) . As we can see, the only important source of background events is the 
SM electroweak processes contributing to the same final state. 

Further enhancement of the signal from the anomalous contributions can be obtained 
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background /cut 


[20 GeV] 






C2H2Z) 


<EU)-(E7J) xp surv 


xp surv 


IRED+- (QCD) 


20.0 


1.12 


0.26 


0.26 


0.058 


0.035 


IRED+- (EW) 


4.4 


0.30 


0.24 


0.24 


0.14 


0.089 


tt 


217. 


6.96 


0.0306 


0.0306 


0.0069 


0.0068 


tij 


1860. 


73.8 


8.88 


0.776 


0.175 


0.158 


tijj 


682. 


77.2 


2.21 


0.0140 


0.0032 


0.0031 


Anomalous o"oo 


2710 


1710 


1310 


1310 


786 


758 



TABLE I: Cross sections in fb to illustrate the effect of cuts for the pp — * jje p^vv production. 
The column marked as (|22 |1 — 1|25 |) [20 GeV] shows the total cross sections after applying out basic 
cuts with a relaxed p™ in = 20 GeV. In computing the SM cross section the contribution from 
a light Higgs boson with = 120 GeV is included. The cross sections (given in fb) do not 
include the forward jet and lepton detection efficiencies and they are obtained for choice CI of the 
renormalization and factorization scales. 

by studying the transverse momentum of the produced leptons as demonstrated in Fig. |3] 
which shows the lepton transverse momentum distribution for the anomalous contribution 
o"oo an d for the SM background. As we can see, the background is peaked toward small 
lepton transverse momenta while the anomalous contributions leads to the production of 
leptons with a higher transverse momentum. Consequently, for processes leading to final 
state leptons with the same charge, we define our signal region by tightening the p l T cut 

Vt > 100 GeV. (31) 



background / cut 


(ZZHZ5) 


CZhCSt + dSU) 


+ <EU>] xP surv 


IRED++ (QCD) 


0.07 


0.004 


0.0009 


IRED++ (EW) 


1.11 


0.105 


0.063 


Anomalous coo 


2250. 


1470. 


880. 


IRED (QCD) 


0.025 


0.001 


0.0002 


IRED (EW) 


0.365 


0.046 


0.028 


Anomalous coo 


536. 


334. 


200. 



TABLE II: Effect of cuts for pp — ► jjepvv production. The cross sections (in fb) do not 
include the forward jet and lepton detection efficiencies and they are obtained for choice CI of the 
renormalization and factorization scales. 

We present in Table ILTT1 our final results for the coefficients Obckj of Eq. (J2U|) after 

cuts <p22|l - (|2Bjl for pp -> jje^i^vv and cuts + (jHU) for pp -> jje ± ^ ± vv. The 
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^0.018 

o 

^0.014 




100 200 300 400 500 

Pt, (GeV) 



FIG. 5: Normalized lepton transverse momentum distribution for IRED++ background (solid 

line), IRED background (dashed line), W + W + a o (dot-dashed line), and W~W~ a 00 (dotted 

line). We assumed mh = 120 GeV and applied cuts 



results include the effect of the veto survival probabilities as well as the forward jet and 
lepton detection efficiencies. They were obtained for choice CI of the renormalization and 
factorization scales. 



D. Estimating the Backgrounds 

As mentioned above, constraining the quartic gauge boson couplings is essentially a count- 
ing experiment and the sensitivity of the search is thus determined by the precision with 
which the background rate in the search region can be predicted. Since the signal selection 
is demanding, including double forward jet tagging and central jet vetoing techniques whose 
acceptance cannot be calculated with sufficient precision in perturbative QCD, the theoret- 
ically predicted background can vary up to a large factor. Though the QCD corrections to 
the irreducible EW processes seem to be modest |22| the same is not guarantee for the QCD 
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scenario 


channel 


Cbck 


O"0 


0"1 


°oo 




O"01 




PP ~ 


-> e ± ^uujj 


0.067 


— 


— 


300 


655 


822 


m h = 120 GeV 


PP ~ 


-> e + n + uujj 


0.029 


-0.46 


-0.20 


400 


94 


380 




PP ~ 


-> e~n~uujj 


0.045 


-0.11 


-0.04 


91 


21 


87 




PP ~ 


-> e ± n^uujj 


0.07 


1.3 


2.1 


300 


655 


822 


No light Higgs boson 


PP ~ 


-> e + n + vvji 


0.017 


-4.9 


-2.3 


400 


94 


380 




PP ~ 


-> e~fi~uujj 


0.017 


-1.2 


-0.54 


91 


21 


87 



TABLE III: Cross sections (in fb) for the different terms in Eq. 1)20(1 . The results include the effect 
of the veto survival probabilities and the forward jet and lepton detection efficiencies. They are 
obtained for choice CI of the renormalization and factorization scales. 

backgrounds. Here the situation is analogous to the Higgs boson production in WBF where 
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the backgrounds must be also estimated from data 

We demonstrate the large QCD uncertainties in Fig. where we plot the value of the ttj 
cross section after cuts (|22|) -(|28| ) for different choices of the factorization and renormalization 
scale. Moreover, we should also keep in mind that the narrow width approximation used 
by us has a discrepancy with respect to the full LO amplitude calculation of the order 
of 10-20% 0]. The obvious conclusion is that in order to obtain a meaningful estimate 
of the sensitivity the background levels need to be determined directly from LHC data. 
Fortunately, a sizable sample of jje ± ix T vv and jje y^vv events will be available if some 
of the cuts are relaxed. In this way the background normalization error can be reduced by 
considering a larger phase space region as a calibration region. The background expected in 
the signal region is then obtained by extrapolation of the measured events in the calibration 
region to the signal region. This procedure introduces also an uncertainty, which we denote 
as QCD-extrapolation uncertainty, due to the extrapolation to the signal region. However, 
as we will show, these uncertainties are smaller than the overall normalization uncertainty. 

Using the results in Fig. 0] we see that we can define the calibration region used to estimate 
the background for jje ± y^vv as the one complying with cuts (}2"2"|) - (f2"%j) and 

M™ w < 800 GeV. (32) 

Equivalently from the results in Fig. EJ we find that one can define the calibration region 
used to estimate the background for jje : /i ,vv as the one within cuts (j22|) -(j2T 



30 < p e T < 100 GeV. (33) 

As a measure of theoretical uncertainty associated with the extrapolation from the cali- 
bration to the signal regions, we study the ratio of the cross sections in the signal region and 
the calibration region as a function of £, the scale factor for the four different renormaliza- 
tion scale choices /zr = ^y° R listed above. In this way we define for the final state exhibiting 
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FIG. 6: Cross section for ttj after cuts l|22 |) -(|28j l as a function £, where //r = ^fJ^ for the four 
choices choices of the factorization and renormalization scale defined in the text. 



opposite charge leptons (jje ^vv 



<?b^{M™ w > 800 GeV) 
a hck (M^ w < 800 GeV) 



where in the evaluation of these ratios we have added the electroweak and QCD contri- 
butions from all background sources taking into account the corresponding veto survival 
probabilities. On the other hand, for the final state jje^/j^uu that exhibits same charge 
leptons we define 

q b ck(^ > 100 GeV) 
ss <x bck (30 < p e T < 100 GeV) ' 1 ' 

where we have added the contributions from both signs. 

We depict in Fig. [7| the £ dependence of R os which shows that the extrapolation uncer- 
tainty is at a tolerable level (~ 15%) being much smaller than the normalization uncertainty. 
The corresponding extrapolation uncertainty for the processes with same sign leptons is 
smaller by a factor of 2 because of the QCD background is small. 

Altogether the total expected uncertainty in the estimated number of background events 
has two sources: the theoretical uncertainty associated to the extrapolations from the cali- 
bration region (<5 D ck,th) and the statistical error associated to the determination of the back- 
ground cross section in the calibration region (5 D ck,stat)- This last one is slightly different for 
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FIG. 7: The ratio R os is shown as a function of £, where hr = f° r the four choices of 
renormalization and factorization scales given in the text. 

the case of light or no light Higgs boson because of the slightly different number of events 
from the SM irreducible background. Assuming an integrated luminosity of 100 fb _1 we find 

£bck,th,os = 15% , £bck,stat,os = 34 ^ ' ^bcdtat.os = 31 % > ( 36 ) 
^bck,th,ss = 7.5% , ^bck.stat.ss = 22 % J ^bcdtat,ss = 21 % > ( 37 ) 

where we denoted by the superscript "lin" ("non-lin") the case with (without) a light Higgs 
boson. In addition to these uncertainties considered here there are also experimental sys- 
tematic uncertainties, which are sizable for the Higgs boson searches j^, however they do 
require a full detector simulation which is beyond the scope of this work. 



VI. RESULTS AND DISCUSSION 

In order to obtain the attainable sensitivity to deviations of the SM predictions of the 
quartic gauge boson couplings we assumed an integrated luminosity of 100 fb -1 and that 
the observed number of events in the different scenarios is compatible with the background 
expectations for the choice CI of the renormalization and factorization scales both in the 
signal (Af data ) and in the calibration (Np data ) regions, i.e. 

-^Vlata = -^i,bck,Cl aI1 d A^data = ^bck,Cl > ( 3 §) 
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where we denote the process pp — *> jje jj^uu by i = os and the sum of pp — > jje+^vv and 
pp — ► jje~fi~vv by i = ss. 

Deviations from the SM prediction for the four gauge boson vertices manifest themselves 
as a difference between the number of observed events and the number of background events 
estimated from the extrapolation of the background measured in the calibration region 

(^?back)> that is > 

^i.data — -^Ys.back > (39) 

where Nf hack = Ri Af c data . Notice that (JH8|) implies that we are assuming that no departure 
of the SM predictions has been observed neither in the control region nor in the signal one. 
The statistical error of the number of anomalous events is 

°f,stat = ^fdata + ( R i ^Pdata^bck.stat) 2 (40) 

where the first term is the statistical error of the measured number of events in the signal 
region and the second term is the error in the determination of the background in the signal 
region due to the statistical error of the background measurement in the calibration region, 
<5bck,stat- The extrapolation uncertainty introduces an additional error 

0-j,th = Ri ^i,data^bck,th • (41) 

Both errors can be assumed to be Gaussian and we combine then in quadrature. 

Given our definition of the signal (|39|. the errors (|40j ) — (j41j ) . and the parametrization of 
the cross section in Eq. (|20j) we can easily obtain the attainable limits on any combination 
of quartic anomalous coefficients. We exhibit in the upper panels of Fig. |H] the 99% CL 
exclusion region in the plane fo versus fi (left) and «4 versus a 5 (right) for each channel % 
independently, and for the combination of both (full region). As we can see the same sign 
processes present a very strong correlation between both couplings while the correlation is 
somewhat smaller for the case of the processes with opposite sign leptons. As a consequence, 
the final allowed regions are rather "compact" and meaningful sensitivity bounds can be 
derived. 

In the lower panels of Fig. |S] we plot the \ 2 as a function of individual couplings, under 
the assumption that only one anomalous parameter is non-vanishing. From this we find 
that for the case with a light Higgs boson of = 120 GeV 

- 22 < -^(TeV" 4 ) < 24 , (42) 
-25 < -£_(TeV" 4 ) < 25 , (43) 

at 99% CL. In models without a light Higgs boson we get the following 99% CL bounds 

- 7.7 x 10~ 3 < ct 4 < 15 x 10" 3 , (44) 
-12 x 10~ 3 < ct 5 < 10 x 10~ 3 . (45) 
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These results represent an improvement of more than one order of magnitude over the 
present sensitivity from indirect effects in low energy observables (|18|) and (|19j) . Notwith- 
standing, for the case in which a light Higgs boson is found and the gauge theory is linearly 
realized, they do not reach the expected natural order of magnitude fi ~ 0{1) for new 
physics scale above 1 TeV. On the other hand, for scenarios without a light Higgs boson a 
natural order of magnitude of the anomalous couplings a, in a fundamental gauge theory 
is g 2 v 2 /A 2 0, since the quartic anomalous interactions can be generated by tree diagrams. 
Thus, we might expect that the size of the a's should be of the order of M§/A 2 ~ 2 x 10 -3 
which is close to the attainable sensitivity that we obtain from our analysis. 

It is also interesting to notice that the achievable sensitivity at the LHC is close to 
the recently derived lower bounds based on the usual analytical properties associated with 
causal, unitary theories 0|. The lack of observation of an anomalous coupling and a§ 
below that bound, would indicate the breakdown of some of these basic properties of the S- 
matrix. In particular, as pointed out in Ref . , since String theory is designed to produce 
S"-matrix with these properties and therefore, the experimental verification of those bounds 
could be use to falsify string theory. 
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APPENDIX A: DIMENSION 8 EFFECTIVE OPERATORS 

We list here the parity conserving effective Lagrangians leading to pure quartic couplings 
between the weak gauge bosons assuming that a Higgs boson has been discovered, that is, 
employing the linear representation for the higher order operators. Denoting by $ the Higgs 
doublet and by U an arbitrary SU(2)l transformation, the basic blocks for constructing the 
effective Lagrangian and their transformations are: 




$ , that transforms as $' = U$> 
, that transforms as D' & = UD^ 

L — > tnat transforms as W = UW^ U U^ 



(Al) 
(A2) 



(A3) 



B, 



that transforms as B' 



B, 



(A4) 
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where W^ v is the SU(2) L field strength and B^ v is the U(l)y one. The covariant derivative 
is given by = (fy - igWfe - ig'B„\)<$>. 

The lowest dimension operator that leads to quartic interactions but does not exhibit 
two or three weak gauge boson vertices is dimension 8. The counting is straight foward: 
when can get a weak boson field either from the covariant derivative of <£> or from the field 
strength tensor. In either case the vector field is accompanied by a VEV or a derivative. 
Therefore genuine quartic vertices are of dimension 8 or higher. 

There are three classes of such operators: 



a. Operators containing just 

The two independent operators in this class are 



£s,o — 


(D^D u <f> 


X 




£>s,i = 




X 


(A,$) f 


b. Operators containing D^<& and field strength 




The operators in this class are: 










£m,o = 


Tr 




X 




£m,i — 


Tr 




X 





£m,2 
£-M,3 



[B^B vff ] x [(Dp^D^ 



x B Pv 
x B^ 



(A5) 
(A6) 



(A7) 
(A8) 
(A9) 
(A10) 
(All) 
(A12) 
(A13) 
(A14) 



c. Operators containing just the field strength tensor 

The following operators containing just the field strength tensor also lead to quartic 
anomalous couplings: 
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= Tr 




x Tr 


w af3 w a ? 


(A15) 




= Tr 




x Tr 




(A16) 




= Tr 




x Tr 


W Pv W va 


(A17) 




= Tr 




X Bj3 v 


(A18) 




= Tr 




x B Pu 


(A19) 


*~i ,b 


= Tr 




x B ap B^ 


(A20) 




= Tr 




x B^B™ 


(A21) 




= Tr 




x B pv B va 


(A22) 




= B^B^ B af3 B 


a/3 






(A23) 


£t,9 








(A24) 
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FIG. 8: Sensitivity bounds on the coefficients of the anomalous quartic gauge boson operators for 
the case with a light Higgs boson (left panels) and with no light Higgs boson (right panels). The 
upper panels show the 99% CL allowed regions for the different channels. The lower panels show 
the dependence of the x 2 function for the different channels assuming that only under that only one 
anomalous parameter is non-vanishing. The dashed (dotted) lines corresponds to pp — ► jje^^vv 
(pp — ► jje^ ijl^vv) while the combined analysis is indicated by the filled region and solid lines. 
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